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Abstract 



Recently Drummond and Hillery [Phys. Rev. A 59, 691 (1999)] presented 
a quantum theory of dispersion based on the analysis of a coupled system 
of the electromagnetic field and atoms in the multipolar QED formulation. 
The theory has led to the explicit mode-expansions for various field-operators 
in a homogeneous medium characterized by an arbitrary number of resonant 
Z^ ' transitions with different frequencies. In this Comment, we drawn attention 

f^ . to a similar multipolar study by Juzeliiinas [Phys. Rev. A 53, 3543 (1996); A 

55, 929 (1997)] on the field quantization in a discrete molecular (or atomic) 
- ^ medium. A comparative analysis of the two approaches is carried out high- 

o ■ lighting both common and distinctive features. 
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Recently Druminond and Hillery [jl[ have presented a quantum theory of dispersive elec- 
tromagnetic modes in a medium characterized by an arbitrary number of resonant transitions 
with different frequencies. The radiation field and the matter was assumed to constitute a 
single dynamical system. In the homogeneous (plane wave) case, explicit expansions have 
been obtained for the field-operators in terms of the operators for creation and annihilation 
of polaritons (elementary excitations of a combined system containing the radiation and the 
matter). Yet, the authors of the paper |Q have overlooked a closely related study |^,^ that 
also provides explicit mode-expansions for various quantized fields, such as the macroscopic 
displacement field, given by (in the Heisenberg picture) [Q: 
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where P^^ ^ (P^ ^ _>^) is the Bose operator for creation (annihilation) of a polariton character- 
ized by a wave- vector k, a polarization index A, and also an extra index m = 1,2, ...,M + 1 
that labels branches of polariton dispersion, M being a number of excitation frequencies 
accommodated by each molecule forming the medium. Here also e^'^^ (k) is a unit polar- 
ization vector, Vo is the quantization volume, cu^ = ck/n^'^^ is the polariton frequency, 
^("i) = fi (uj^'] is the refractive index (calculated at (^^ ), and f^™'^ = dcuj!^' /dk is the 
branch- dependent group velocity. 

In what follows we compare the formalism by Drummond and Hillery IQ to that by 
Juzeliiinas p|-§[ highlighting common and distinctive features of the two approaches. Both 
studies consider a similar coupled system of the radiation field and the matter, exploiting the 
same multipolar formulation of Quantum Electrodynamics (QED). (Yet, different techniques 
have been employed to represent the operators of interest in terms of the normal polariton 
modes.) Moreover, in either approach an arbitrary number of transition frequencies (of 
electronic or vibrational origin) has been included for the material medium. As a result, 
the above mode-expansion (|l])-(||) reproduces the same functional dependence on the group 
velocity and the refractive index as the corresponding expansion for the macroscopic dis- 
placement field in the one-dimensional case, given by Eq.(5.17) of ref. |l|]. Such a functional 
dependence is also consistent with earlier narrow-band Langrangian approach by Drum- 
mond 1^. Furthermore, the mode-expansion (|l|)-(@) is equivalent to the three-dimensional 
displacement operator derived by Drummond and Hillery [|1],0, as long as the spatial dis- 
persion is neglected in the corresponding Eq. (8.24) of ref. ||l|]. The same holds for other 
field-operators, such as the operator for the transverse electric field whose mode-components 
are related to Eq.(^ via a relationship of the classical type gj: d,^,^ (r) = €oel"^^e^^ (r), the 

emerging relative dielectric permittivity e^™^ = m^^M being a branch-dependent quantity. 
This is also in agreement with the ref. |l[]. 

Both studies take special care in making sure that the (equal-time) commutation rela- 
tionships preserve between various operators in their diagonal (polariton) representation. 
The study [^^ has checked the commutation relationships involving the field-operators d , 



e-*-, a^, p-*- and h . Their correctness appears to be ensured by the following equalities 
iii: 

E t;("^)n(™) = c and E w^'^V^^"'^ = c (3) 
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Drunimond and Hillery |I| have also exploited one of the these equalities in analyzing the 
commutation relationships. In addition, a few more complex equalities have been presented 
when dealing with the material operators. 

A distinctive feature of the approach by Drummond and Hillery [Q is inclusion of spatial 
dispersion, i.e. dependence of the dielectric permittivity not only on the frequency, but also 
the wave- vector [0. This does not alter the form of the above relationships, as well as the 
mode-expansions of other macroscopic field-operators presented earlier p|, yet the meaning 
of group velocity is to be modified [^ . It is noteworthy that the spatial dispersion has been 
included in the 'effective mass' approximation through an extra differential term featured 
in Eqs.(6.2) and (6.8) of ref. |^: The term represents the coupling between infinitely close 
dipole-oscillators comprising the continuous dielectric medium. One might argue that such 
an approach is not fully consistent with the spirit of multipolar QED formulation in which 
there is no direct coupling between the dipole-oscillators |Q. However, this is perhaps the 
only way to include the spatial dispersion within the continuous model of the dielectric 
considered. 

On the other hand, the approach p|-^ assumes the matter to be discrete, the constituent 
molecules (dipole-oscillators) forming a cubic lattice. Here the effects of spatial dispersion (as 
well as other intermolecular coupling) are contained in the initial multipolar Hamiltonian for 
the radiation field coupled to the discrete medium, all interatomic coupling being mediated 
exclusively via the transverse virtual photons P]. As a result, the spatial dispersion is 
implicit in the general analysis of ref. [Q] up to Eq.(3.40). The effect has been omitted in 
the subsequent long wave-length approximation made in Eq.(4.1) of ref. 0] followed by the 
explicit results. In fact, the spatial dispersion does play a minor role for the optical (photon- 
like) modes characterized by small wave-vectors. One can recover the spatial dispersion in 
a relatively straightforward manner using the discrete model 0, however this is out of the 
scope of the present Comment. Note that in contrast to the continuous approach |]I[, the 
spatial dispersion is characterized exclusively by the microscopic parameters of a discrete 
system, there being no need to introduce an extra parameter describing the effect. 

Using a discrete approach, one can also recover the local field effects from first princi- 
ples. In doing this, the theory @-Q treats systematically the Umklapp processes playing 
an important role in the multipolar QED formulation. As a result, the required local- 
field corrections emerge intrinsically in the refractive index and the group velocity entering 
the mode-expansions of the field-operators [0-^. Furthermore, the discrete approach al- 
lows us to consider not only the macroscopic operators 0, but also the operators for the 
local 10-^ and microscopic P fields. For instance, the mode-components of the local dis- 
placement operator are related to those for the macroscopic displacement operator as 0: 

^i.,m i^d = v^r'^j (^r""^ + 2) /3 dj^,„ (r^). This appears to be very helpful for the analy- 
sis of various molecular-radiation processes in dielectric media 0,^ , such as the spontaneous 
emission. 

Finally, Drummond and Hillery 0] have pointed out that the solution to the eigenvalue 
equation tu^ = c^k"^ /n^ given by Eq.(7.5), 'is unique for any given modal frequency, but has 



forbidden regions which indicate a resonance, or absorption, band'. This would be absolutely 
true for the spatially non- dispersive media, as illustrated in fig.2 of ref. |0]. However Eq.(7.5) 
of ref. [^ contains the spatial dispersion (in the 'effective mass' approximation), so that 
n = n {u, k). Inclusion of such a spatial dispersion yields more than one value of k for certain 
frequencies P , the additional solutions representing the exciton-like modes characterized by 
much larger k. 
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